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We describe a wave-mechanical implementation of classically chaotic n-disk scattering based on 
thin 2-D microwave cavities. Two, three, and four-disk scattering are investigated in detail. The 
experiments, which are able to probe the stationary Green's function of the system, yield both fre- 
quencies and widths of the low-lying quantum resonances. The observed spectra are found to be in 
good agreement with calculations based on semiclassical periodic orbit theory. Wave-vector auto- 
correlation functions are analyzed for various scattering geometries, the small wave-vector behavior 
allowing one to extract the escape rate from the quantum repeller. Quantitative agreement is found 
with the value predicted from classical scattering theory. For intermediate energies, non-universal 
oscillations are detected in the autocorrelation function, reflecting the presence of periodic orbits. 
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I. INTRODUCTION 

The system of n-disks on a plane is perhaps the sim- 
plest and paradigmatic example of a chaotic scattering 
system The classical differential cross-section of 

the system with n > 2 disks on a plane is singular, with 
singular points forming a Cantor set jl],[|. Thus, the 
investigation of chaotic scattering geometries naturally 
brings in the capability to access the intrinsic fractal na- 
ture of the underlying classical repeller. 

From the perspective of quantum-classical correspon- 
dence, the wave-mechanical counterpart of this "fractal 
pinball game" is even more intriguing. Semiclassical at- 
tempts to treat this problem led to important theoretical 
advances and to the development of sophisticated semi- 
classical tools, notably cycle expansion, that represents 
one of the more productive applications of periodic orbit 
theory j5|. 

From a broader perspective, the n-disk scattering can 
be regarded as the prototype of an open quantum chaos 
system. For closed quantum chaotic systems (billiards 
in particular), theoretical and experimental results are 
available on spectral eigenvalues (which are purely real), 
and a good understanding of both universal and non- 
universal features has been reached ||. Open chaotic 
systems are a special class of open quantum systems, 
whereby bound states of a closed geometry are converted 
to long-lived metastable states due to coupling to con- 
tinua. Accordingly, eigenvalues are intrinsically complex 
and their universal behavior is currently a subject of con- 
siderable interest 0. 

Investigation of n-disk scattering was originally stimu- 
lated by the attempt of modelling unimolccular reactions 
j|. However, the system turns out to be a good exem- 
plification for a variety of physical situations, ranging 
from crossroad geometries for semiconductor devices, to 
electromagnetic and acoustic scattering, and heavy-ion 



nuclear reactions || . 

In spite of extensive theoretical analysis, there have 
been few direct experimental implementations of such 
chaotic model geometries so far. In this paper we present 
the results of a systematic experimental study of the 
quantum resonances and decay characteristic of 2-D scat- 
tering repellers. The experiments utilize thin microwave 
geometries where, under appropriate conditions, the un- 
derlying Helmholtz equation maps exactly into the time- 
independent Schrodinger equation in two dimensions. 
Experiments were carried out for n = 1,2,3,4, as well 
as for large n = 20, the latter corresponding to the so- 
called random Lorentz scatterer. A concise account of 
the main results has been presented in H. 

The experiments, which are able to access the station- 
ary Green's function of the system, yield the frequencies 
and the widths of the low-lying quantum resonances. We 
carried out semiclassical calculations of such resonances, 
which are found to reproduce the measured spectra rea- 
sonably well. Our experiments enable us to explore the 
role of symmetry in a unique way by probing different 
irreducible representations of the symmetry group of the 
scatterer. The experimental data are used to identify 
both universal and non-universal signatures of the classi- 
cal chaos in the transmission spectra, through measures 
such as the spectral (wave- vector) autocorrelation func- 
tion. 

Correlation functions are a valuable tool to extract key 
information on the spectral properties of the system. In 
mesoscopic conductors, measurements of the magnetic 
field correlation of the conductivity have provided unique 
insight into the manifestations of the chaotic classical dy- 
namics at the quantum level [Tcfl . Correlations of wave 
functions in chaotic systems have been considered re- 
cently as a probe for classical ergodicity predictions [^l) . 
In the present experiments, we take advantage of our 
ability to vary the wave-vector of the system to directly 



1 



access energy correlations, that are difficult to extract in 
semiconductor microstructures. The small k (long time) 
behavior of the resulting autocorrelation provides a mea- 
sure of the quantum escape rate, which is shown to be in 
good agreement with the corresponding classical value. 
For large k (short time), the contribution of periodic or- 
bits is evidenced through non-universal oscillations of the 
autocorrelation function. 

The content of the paper is organized as follows. Af- 
ter describing the experimental realization in Sec. II, 
the connections between electromagnetic and quantum 
mechanical scattering, and between the observed trans- 
mission amplitude and the stationary Green's function 
are elucidated in Sec. III. In Sec. IV, the essential back- 
ground on semiclassical theory is recalled, with empha- 
sis on the role played by symmetry properties. The re- 
maining sections are devoted to the presentation and the 
interpretation of the experimental results: in Sec. V, 
we present a comparison between the measured traces 
and the semiclassical predictions for various scattering 
geometries probed in the experiment, while in Sec. VI we 
focus on the analysis of spectral autocorrelations, from 
which the wave-mechanical escape rate can be extracted. 
An overall discussion of the results and their implications, 
together with an outlook on open issues and future work 
are presented in the conclusive section. 

II. EXPERIMENTAL REALIZATION 



edge of E z alone. The output signal can be picked up 
by both directly coupling to the axial electric field E z 
(electric coupling) or by measuring the associated trans- 
verse magnetic field H (magnetic coupling). In the first 
configuration, a microwave pin is inserted through appro- 
priate holes drilled on the top plate of the cavity, while 
magnetic coupling is established by suitably positioning 
a second loop perpendicular to the plates. Note that, for 
a sufficiently small loop area, the measured signal is sim- 
ply proportional to the electric field value at the pick-up 
probe location in both cases. Data from both electric and 
magnetic coupling configurations are used in the subse- 
quent analysis. 

The resonances were probed in the transmission mode 
by sweeping the frequency in the range between and 
20 GHz. All the measurements are carried out using an 
HP8510B vector network analyzer measuring the com- 
plex transmission coefficient 5*21 (/) = (X + iY)(f) of the 
coax plus scatterer system. It is crucial to ensure that 
there is no spurious background scattering due to the fi- 
nite size of the system. This was verified carefully as well 
as that the effects of the coupling probes were minimal 
and did not affect the results. Also, we stress that dissi- 
pation effects due to finite conductivity of the walls are 
entirely negligible in the present experiments. 

III. GREEN'S FUNCTION AND 
EXPERIMENTAL S-PARAMETERS 



The experiments were carried out in thin microwave 
structures consisting of two highly conducting Cu plates 
spaced d ~ 6 mm and about 55 x 55 cm 2 in area. A pic- 
ture of a typical experimental setup is displayed in Fig. 1. 
Similar experiments in closed geometries have provided 
direct observation of scars [O, enabled precise tests of 
eigenvalue statistics fl3| , and allowed experimental stud- 
ies of localization effects [[l4| . Disks and bars also made 
of Cu and of thickness d were placed between the plates 
and in contact with them. Disks of radius a = 2 and 
a = 5 cm were used in the experiments. In order to simu- 
late an infinite system, microwave absorber material was 
sandwiched between the plates at the edges. Microwaves 
were fed into the system by inserting a loop-terminating 
coaxial line in the vicinity of the scatterers. Since we 
ensure that the relation / < f c = c/2d = 25 GHz is 
valid for all operating frequencies (here k = 2irf /c, c be- 
ing the speed of light), the only allowed electromagnetic 
modes are a class of TM (Transverse Magnetic) modes 
with no variation along the direction z orthogonal to 
the plates. The nonvanishing field components for these 
states are the axial component of the electric field E z , 
and transverse components of the magnetic field H x , H y , 
or H = (l/ik)z x WE Z . Thus, essentially, the geome- 
try of the experiment can be regarded as 2-dimensional, 
the whole field configuration being accessible from knowl- 



A. Electromagnetic vs. quantum mechanical 
scattering 

According to Maxwell equations, the problem of sta- 
tionary scattering of electromagnetic waves from perfect 
metallic conductors in a 2-D plane is defined mathemat- 
ically by the Hclmholtz equation for the field component 

(y 2 + k 2 )E z (r) =0, (1) 

for f = (x,y) outside the scattering region (i.e., between 
the disks), supplemented by Dirichlet boundary condi- 
tions on the perimeter of any scatterer: 

E z {f) = 0, r £ perimeters, (2) 

implying, for any geometry, a vanishing electric field in- 
side the disks. The Helmholtz equation is formally iden- 
tical to the time-independent Schrodinger equation with 
^ = E z . Since, in addition, the same boundary condi- 
tions apply to both the axial electric field and the wave- 
function, the problem is equivalent to quantum mechan- 
ical scattering from hard disks. This correspondence is 
exact in the whole range of frequencies below the cutoff 
frequency f c that are considered in the experiment. 

According to quantum scattering theory, the station- 
ary properties of the scattering process are characterized 
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completely by the collection of S-matrix elements, ex- 
pressing the transition amplitudes between asymptotic 
incoming and outgoing states HI01- ^he ^-matrix ex- 
hibits poles as a function of the wave vector (or en- 
ergy) of the incoming wave. These poles, corresponding 
to scattering resonances, occur in the complex fc-plane, 
k = k' + ik" . The imaginary part is interpreted as a 
resonance width, implying that scattering resonances are 
associated with metastable states having a finite lifetime. 
Note that, in contract to closed-cavity geometries where 
non-zero widths only arise due to the finite wall conduc- 
tivity [jl6|, widths have here a purely geometric origin, 
the openness of the system preventing a stable trapping 
of the states within the scattering geometry for arbitrar- 
ily long times. We shall use the convention k = s — ij 
for a given resonance, i.e., s — k' to denote the position 
of the resonance along the fc-axis and 7 = — k", k" < 
to denote half-widths in the wave-vector domain. 7 is 
related to the observed resonance width A/ in the fre- 
quency spectrum by A/ = (c/2n) Ak — (c/2n) 2j. 

The correspondence, pointed out above, between the 
equations of motion and the boundary conditions sat- 
isfied by electromagnetic (em) and quantum-mechanical 
(qm) stationary scattering allows one to establish a di- 
rect mapping between the S'-matrix spectral properties 
in the domain of complex wave vectors: 

(k' + ik") EM = {k' + ik") QM . (3) 

This mapping enables us to study the quantum proper- 
ties of the n-disk system. Note that quantum-mechanical 
(half) widths T/2 for unstable states are typically defined 
in the energy domain as E = E' + iE" = e - iT/2, E" < 
0, the relation between (complex) pseudo-energies and 
(complex) wave-vectors being E = (hk) 2 /2m as in the 
real case |L7|, @. Thus, E" = {fi 2 /m)k'k". How- 
ever, due to the difference between the non-stationary 
propagators for em and QM problems, care should be 
taken when introducing and comparing associated quan- 
tities like decay rates or lifetimes which are defined in 
the time domain. By exploiting the time-dependent 
Schrodinger equation, quantum decay rates for square 
amplitudes ^(i)! 2 are found as inverse lifetimes 1/r = 
T/h= \2E"\/h. One gets r -Y/h= 1, expressing Heisen- 
berg's uncertainty principle. From the full wave equation 
for a monochromatic field, the decay rate of the field am- 
plitude \E z (t)\ 2 is given by 1/r = 27c = c\2k"\. Accord- 
ingly, t • 27c = 1. In spite of such a difference, the wave 
mechanical spectrum will allow us to extract a quantity 
(the classical escape rate) which is related to the aver- 
age trapping time of the underlying classical trajectories 
inside the scattering region. 



B. Green's function 

We turn now to outline the relationship between the 
measured response S21 and the two-point Green function 
of the chaotic billiard. The equation of motion satisfied 
by the relevant field component E z in the presence of a 
monochromatic driving at frequency u) is 

(V a -^)j5.(*',*)=^(*0e-*-*, (4) 

r denoting 2-D coordinates in the plane and Eq (r) being 
the associated amplitude for waves fed into the cavity 
with wave vector k = lo/c. For a stationary field config- 
uration, the above equation reads 

(W 2 + k 2 )E z (r) =E (r). (5) 

By definition, the so-called Green's function (or propa- 
gator) associated with the wave equation differential op- 
erator satisfies the equation 

( y2 ~ G ^ r ; t, t ) = S(r- r )5(t - t ) , (6) 

for fixed ro and to- The Green's function G(r, ro) t, to) is 
identical to the wavefunction generated by the system at 
point f and time t in response to a delta-like excitation 
applied at point r*o and time to. For a time-homogeneous 
system, G(r, ro; t, to) = G(r, ro; t — to), and Fourier trans- 
formation of Eq. (^) gives 

(V 2 + k 2 ) G(f, r ; k) = 6(f- f ) , (7) 

where the stationary Green's function has been intro- 
duced: 

G(f, r ; t - io ) = i- / G(f, r Q ; k) e **(*-*<0 d k . (8) 

Knowledge of the Green's function allows one to calculate 
the response of the system to a generic excitation Eq (r) 
as a convolution integral E^r^) — J v G(r 2 , f; k)Eo{r) dr. 
In particular, if a point-like probe is placed at f\ to ex- 
cite the system, E {r) = Eq(ti) 8{r — r\), we simply get 
E z {r 2 ) = G(r 2 ,r 1 ;k)E Q {r 1 ). 

The connection with the observable quantity of our 
experiments is established by recalling that the trans- 
mission coefficient S21 = V ou t/Vi n , where Vi n and V ou t 
are linear responses to the electric field at the probe lo- 
cations fi (input) and r 2 (output) respectively. Note 
that S21 can be related to the S'-matrix element de- 
scribing transmission between waves entering through 
antenna 1 and going out through antenna 2 Wc 
have Vout = ctE z (r 2 ), V m = f3E {ri), a and (3 denoting 
impedance factors that are characteristic of the coax lines 
and the analyzer, and are generally slowly- varying func- 
tions of frequency. Thus we get S21 = {&■/($) G(r*2, rf; k), 
which we can write in the form 
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S 21 (k)=A(k) G(r 2 ,r i; fc). 



(9) 



In a formal analogy with the closed-system case, the 
Green's function can be expressed in terms of a gener- 
alized eigenfunction expansion |l7f| , 



G(r 2 ,n; k) = ^ — 2 p 



(10) 



* n (r) denoting the n-th eigenfunction of an effective 
open-system Hamiltonian, with corresponding adjoint 
eigenstate 3> n (r) and associated (complex) wave vector 
k n . It is worth stressing that the measured Green's func- 
tion is identical with the actual Green's function of the 
open system provided the perturbation introduced by 
the probes is negligible. If this is not the case, a modi- 
fied Green's function is generally obtained, with shifted 
complex poles in the expansion ( |l0| ) (i.e., resonances are 
shifted and broadened due to the presence of the probes) . 
In the experiment, we ensured that the antennas were 
only weakly coupled to the system, with associated per- 
turbation effects less than w 10~ 3 of the observed fre- 
quencies and widths. Under these conditions, the relation 
established above between S21 and the Green's function 
generalizes the derivation for closed microwave billiards 
presented in pj|. 



C. Comparison with mesoscopic conductivity 
measurements 

The connection between so-called transmission func- 
tion, S-matrix, and Green's function is well known in 
the field of mesoscopic transport |l7j]. In particular, we 
recall that the transmission function of a coherent con- 



ductor between two leads 1, 2 is defined as 20 



T / 1 ; in a 

mel,ne2 



(11) 



where the transmission coefficient t mn , which character- 
izes the transmission amplitude between mode m in lead 
1 and mode n in lead 2, is given by 

t mn = ~ifi\/vmv n J dy'dycj>^ n (y')cj) n (y)G(y',x 2 ,y,x 1 ;k). 

(12) 

Here, v m (v n ) and </> m (0n) are the longitudinal velocity 
and the transverse wave function for the mode m in lead 
1 (n in lead 2) respectively, while ^1(2) denote the lon- 
gitudinal coordinate of antennas 1,2. The transmission 
coefficients, which are clearly energy (frequency) depen- 
dent, are directly related to the correspondent elements 
of the S-matrix. In fact, the above expression for t mn is a 
two-probe version of the general connection known as the 
Fisher-Lee relation between the S-matrix and the Green's 



function Q|. The transmission T is related directly to 
the conductance a of the conductor by the Landauer for- 
mula, a — (e 2 /h) T, e denoting the electron charge [Ti| . 

In the case of a point-like excitation and a point-like 
probing of the system, the input / output leads act as zero- 
dimensional tunneling point contacts. Thus, the trans- 
verse dimension of the leads can be neglected and the 
wave functions <j> m (y) m Eq. ( |l2] ) are proportional to 
^-functions, </> m (y') = S(y' - y 2 ), <j> n (y) = <% - Vi)- 
By combining the two expressions above, and putting 
?l = [xx,yi), r 2 = (x 2 ,y 2 ), we get T(k) cx \G(r 2 , fi, k)\ 2 . 
Since the Green's function is directly related, as shown in 
the previous section, to the experimental trace, one can 
interpret 



T(fc)oc|S 2 i(fc)| 5 



(13) 



as a measurement of the two-point conductance. As a 
consequence, direct comparison is possible with theo- 
ries originally developed for electronic micro-structures 
|Io| , pO 22 2^]. Note that in the language of mesoscopic 
conductors, the weak-coupling assumption between the 
system and the point-like leads implies the possibility of 
neglecting, as noticed above, any perturbation effect as- 
sociated with the so-called lead self-energy |l7j . From a 
general point of view, microwave-based implementations 
offer, compared to their solid-state counterparts, the ad- 
vantage of a simple and practically unlimited manipula- 
tion of geometrical properties, along with the possibility 
of changing the wave-vector (and thereby the energy) of 
the incoming waves at will. Lastly, since electromagnetic 
waves are intrinsically non-interacting, a direct analogy 
with the ideal limit of a non-interacting electron gas 
applies and microwave experiments automatically probe 
quantum transport in the ballistic regime. 



IV. SEMICLASSICAL THEORY 

The starting point for the semiclassical derivation is 
the concept of generalized density of states D(E) which, 
according to Balian and Bloch |24[] , is defined as a suit- 
able difference between the density of states of the free 
and of the scattering system. The relative density of 
states is related to the S'-matrix in the following way: 



D ^ = ^r 



fit 



dS(E) 
dE 



(14) 



Therefore, the density D(E) and the S'-matrix share the 
same complex poles. D(E) itself is written as 
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D{E) = lmg{E) 



(15) 



where the function g(E), which is related to the so-called 
quantum Selberg zeta function, is the trace of the differ- 
ence between the Green's function in the presence and in 
the absence of the disks respectively. 
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In the semiclassical limit, g(E) is expressed as a sum 
over periodic orbits i.e., by Gutzwiller's trace formula 
p5[ . The trace function can then be written in terms of 
the Ruelle (-function as [[|, 



d 



9(E) =9o-J2d^ ln C(i/2)+j(-«fc) 



(16) 



where go, which is independent on the energy E, is 
given by the difference between the so-called Thomas- 
Fermi state densities with and without the disks, e.g., 
go = —Aa 2 m/(2h 2 ) in the 4-disk configuration. In Eq. 
(|l6|), the Ruelle (^-function with j running from to oo 
is given explicitly as 



c ( i/2 )+J -Hfc)=n i 1 -(-^ 



/4 



l/2)+J 



(17) 



where k is the wave vector, l p is the length of the periodic 
orbit p , L p the number of collisions of the periodic orbit 
with the disks, and A p the eigenvalue of the so-called sta- 
bility (or monodromy matrix). According to the above 
expressions, the task of finding the scattering resonances 
of the system is reduced to the problem of identifying 
the complex poles of the Ruelle ^-function. Since the 
Ruelle £ -function with j — provides the leading contri- 
butions to the resonances with longer lifetimes (sharper 
resonances), we will only restrict ourselves to this case in 
our subsequent discussion. 

The Ruelle ^-function can be conveniently expressed 
in terms of a so-called Euler product representation, 



(18) 



t p denoting the semiclassical weight for cycle p pq] . In 
practice, the periodic orbit formula for the Ruelle func- 
tion is evaluated by performing a cycle expansion and 
investigating the zeros and radii of convergence as func- 
tions of truncations to cycles of a given maximum topo- 
logical length. In the n-disk problem, if the disks are 
positioned sufficiently far from each other along a ring, it 
is possible to travel between any three successive of them 
and the trapped trajectories can be put in one-to-one 
correspondence with the bi-infinite sequence of symbols 
ujk taken in the alphabet {1, 2,3, ... ,n} with u>k+i ^ &k- 
The latter constraint implies that the topological entropy 
per bounce is equal to ln(n — 1). Thus, the system sym- 
bolic dynamics has a finite grammar and the Euler prod- 
uct (jl^) can be rewritten by separating out a dominant 
fundamental contribution and the remaining corrective 
terms: 



na-*iO=i-£*/-£' 



(19) 



Here, the number of fundamental terms tf is equal to 
the number of symbols in the unconstrained symbolic 



dynamics, while the so-called curvature corrections c r 
represent contributions due to the nonuniformity of the 
system p6|. 

In our experiments, n identical disks are placed on a 
ring with equal space between the nearest neighbors. Ac- 
cordingly, the system is characterized by symmetry point 
group G = C 2v , C 3v , Ci v , C' 5v , C ev ,... for n = 2,3,4, 
etc., and symmetry properties can be used to classify 
the corresponding scattering resonances. In practice, the 
presence of symmetries can also be exploited to simplify 
the cycle expansion and improve its convergence. The 
key concept is to remove symmetry-induced degenera- 
cies between cycles by reducing the dynamics to the so- 
called fundamental domain (FD). The latter is a region 
obtained by ideally replacing the symmetry axes with 
perfectly reflecting mirrors. Global periodic orbits of the 
full system can be described completely by folding irre- 
ducible segments into periodic orbits in the fundamental 
domain. Correspondingly, the Ruelle £-function can be 
factorized in the product over different irreducible repre- 
sentations (irreps) of the symmetry group |26f| . By taking 
into account the degeneracies of the periodic orbits and 
working in the fundamental domain, the Euler product 
(p8|) takes the form 



n( i -^)=ri( i 



t h py 



(20) 



where a modified semiclassical weight has been intro- 
duced: 



In Eqs. (E(H)-(E1 



(21) 



g/hp, g and h p denoting re- 



spectively the order of the full symmetry group G and 
of the maximal subgroup leaving p invariant, p is the 
irreducible segment of the orbit p that corresponds to a 
fundamental domain orbit, Lp is the number of collisions 

with the disk in the fundamental domain and Ap = Ap 
the eigenvalue of the stability matrix in the fundamental 
domain. The calculation of Ap can be easily performed 
by diagonalizing the stability matrix Jp which, for hard- 
disk billiards, reads as |26lE!7| 



.I)-, n 



fc=i 



i h 

1 



1 

Tk 1 



(22) 



Here, for simplicity we rescaled the constant velocity to 
unit value. Ik denotes the length of the fc-th free-flight 
segment of the cycle p, while Tk = 2/a cos (f>k is the de- 
focusing due to the fc-th reflection, occurring at an inci- 
dence angle 4>k- 

By construction, only one irreducible representation 
appears whenever the system is probed in the funda- 
mental domain. Note that, in practice, one can easily 
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switch from one symmetry to another by varying the 
angle 8 — n/n. We recall now some relevant formulas 
for the symmetry configurations examined in the exper- 
iment. 



A. C2 factorization 

The group G = C% v is the appropriate symmetry group 
for two-disk scattering which, without being chaotic, of- 
fers the most important example of a nontrivial inte- 
grable scattering problem. An exact analytical solution 
is available for the classical dynamics |^8[ , while a full 
quantum mechanical calculation of the scattering reso- 
nances has been performed in [^9|,^l). From the stand- 
point of the semiclassical analysis, the symmetry group 
of the periodic orbits is G = C2 , characterizing the trans- 
formation properties of the orbits under the exchange of 
the two disks |p9| . The fundamental domain is the half- 
space containing a single disk. There are just two group 
elements, C2 — {e, P}, e being the identity and P the 
parity operation, and only one periodic orbit. Cycles 
classify according to the two irreducible representations 
A\ (symmetric) and A2 (antisymmetric). Let a and R de- 
note the radius and the center-center disk separation re- 
spectively, with the ratio a = R/a. We have 



C A l = 1 - to 
with to — — exp[i/c(i? - 

A = (<t — 1) + a/cKo 7 



CaI = 1 + to : 

2a)} /VA, and 



2), 



(23) 



(24) 



denoting the eigenvalue of the monodromy matrix indi- 
cated above. 

We thus get the semiclassical scattering resonances as 



(2n — m) 7r 



■In A 



R-2a 



n= 1,2,. 



(25) 



where m = 1 for the j4i-irrep, m — for the ^-hrep. 
In the fundamental domain, only the antisymmetric A2 
representation contributes. 



B. C% v factorization 

The symmetric three-disk pinball is invariant under the 
transformations of the group C^ v . In addition to the iden- 
tity transformation, they include two rotations through 
2tt/3 and 4-7r/3 about the main axis, and three mirror 
reflections around the symmetry axes. The fundamental 
domain is bounded by a disk and the two adjacent sec- 
tions of the symmetry axes acting as mirrors ( 1 / 6th of 
the full space, see Fig. 3, inset). Cs v has two 1-D irreps 
A± and A2 (symmetric and antisymmetric under reflec- 
tions respectively) , and one 2-D irrep of mixed symmetry 



labelled E. The 3-disk dynamical zeta function factor- 
izes into: £ = Ca 1 Ca 2 (e' the contributions of each given 
irreducible representation being given by the following 
curvature expansion j26|: 



Cai — 1 — *o — t± — (toi — toh) 

— [(tool — totoi) + (^011 — tltoi)] 



for the Aisubspace, 

Ql = l + to-ti + (*oi-*o*i) 
— [(tool — totai) — (ton — 



tihi)} 



(26) 



(27) 



for the antisymmetric A2 subspace, and, for the mixed- 
symmetry subspace E, 



Ce 1 = 1 + ti - (4 - t\) + (tooi - Mo) 

+ [to011 + (tool — t\t>o)t\ — tgj + 



(28) 



The representation in the fundamental domain is A2. 
A detailed comparison between the semiclassical predic- 
tions and the exact quantum resonances for the 3-disk 
scattering problem has been reported in @J30]. Exact 
quantum calculation was also performed in 



C. Civ factorization 

The scattering problem of four equal disks placed on 
the vertices of a square is characterized by C^ v symmetry. 
This is a group consisting of the identity, two reflections 
across the coordinate axes, two diagonal reflections, and 
three rotations by angles ir/2,ir and 37r/2. The funda- 
mental domain is a sector delimited by a disk, a portion of 
the corresponding diagonal axis and a portion of the con- 
current coordinate axis ( i.e., l/8th of the full space, see 
Fig. 4b, inset). C^ v has four 1-D irreps, either symmetric 
(Ai) or antisymmetric (A2) under both types of reflec- 
tions, or symmetric under one and antisymmetric under 
the other (Bi, B2), and one 2-D representation E. The £- 
function is factorized as C = Ca 1 (a 2 Cb 1 Cb 2 (e' wnere the 
contributions for the various invariant subspaces have the 
following curvature expansion p6| : 
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+(iooi — totoi) — (too2 — £0^02) — (ton ~ titoi) 
-<023 + *ia2-(*na-ti*i2) (29) 
+ (^012 + £021 ^ ^0^12 — £1*02) • ■ ■ j 
C, B l = l+to+ti-t 2 

-(*oi - *o*i) + (toa - *o*a) + (*ia - *i*a) 

+ (iooi — io<oi) — (^002 — £0*02) + (ion — Woi) 

+ (^022 — t 2 to 2 ) — (*112 — ^1*12) + (*122 — ^2^12) 

— (ioi2 + tei + £0*1*2 — toti2 — t±to2 — *2*oi) • ■ • j 

Cfi 1 = 1 + *2 ~ (*o ~ ^l) + (2*002 — *2*o — 2tn2 — *2*l) 
+ (2*0011 — 2to022 + 2*2*002 — *01 — *02 
+2*1122 — 2*2*112 + *12 — *0*l) ■ ■ ■ ■ 

The representation in the fundamental domain is £?2- 
A detailed comparison between the semiclassical predic- 
tions and the exact quantum resonances for the 4-disk 
scattering problem has been performed by f27fl , 

V. COMPARISON WITH EXPERIMENTAL 
RESONANCES 

The experimental transmission function |S , 2i(fc)| 2 = 
X 2 (k) + Y 2 (k) (Sec. II) can be expressed as a super- 
position of Lorentzian peaks, 

where, as above, Sj, ji denote respectively the position 
and the half-width of the resonances in the fc-domain. 
The parameters c$ are coupling coefficients that depend 
on the location of the two probes and reflect the cou- 
pling between the pick-up antenna with the -E z -pattern of 
a given resonant mode. Semiclassical calculations using 
the appropriate cycle expansion described in the previous 
section were performed for different geometries, leading 
to the real and imaginary part Si and ji of the resonances. 
In comparing with the observed traces, the parameters 
q were set manually to fit the data. For a given scatter- 
ing geometry, the trace S21 used is the average of several 
traces collected at different probe locations in order to 
avoid missing resonances due to the accidental coinci- 
dence of either probe with a node of the wavefunction. 
In general, good agreement is found for the resonant fre- 
quencies of relatively sharp resonances (with typical qual- 
ity factors in the range Q = f / A/ < 50). Broader reso- 
nances with larger imaginary parts are instead not easy 
to distinguish, although all resonances are always con- 
tributing to the transmission function. We examine now 
specific configurations. 

2-disk: For the two-disk scattering, preliminary mea- 
surements were reported in p2[ . We carried out exper- 
iments in both the full and in the half-space geometry, 
with a = 5 cm and R = 40 cm. According to the dis- 
cussion in Sec. IVa, the trace is expected to exhibit 



resonance peaks at regularly spaced locations, f n = n 
GHz, n = 1,2, ... , with a constant width approximately 
equal to Af n ~ 0.29 GHz (from Eq. (§5|)). A typical 
experimental trace is shown in Fig. 2, where we focussed 
on the A2 -resonances between and 20 GHz. The cor- 
responding calculated trace is depicted as a dashed line. 
The agreement is found to be quite good for both the 
resonances and their width. The regularity of such a 
spectrum will manifest clearly in the corresponding au- 
tocorrelation function. 

3- disk: For the 3-disk geometry, we recall that a first 
demonstration of classical chaotic scattering via scatter- 
ing of laser light was presented in f33|j . A typical mi- 
crowave trace for a 3-disk scatterer with a = 4-^/3 is 
presented in Fig. 3. Again, we focus on the fundamental 
domain representation of the scattering geometry, cor- 
responding to resonances with ^-symmetry. The semi- 
classical calculations, which are shown as a dashed line, 
are carried out by using the cycle expansion ( |l7|) with 
8 periodic orbits up to period 4. We verified that they 
accurately reproduce previous calculations on the same 
system [ pM|j3C[ | . For this scattering geometry, comparison 
with the exact quantum mechanical calculations is also 
available |jT| , implying a stringent test for the validity 
of the semiclassical method. According to Fig. 3, the 
overall agreement is qualitatively good, especially for the 
locations of the sharper resonances. 

4- disk: The traces of a 4-disk scatterer with a = 4 in 
the full space and the fundamental domain are shown in 
Fig. 4 top and bottom, respectively. Semiclassical calcu- 
lations (dashed line) were performed by including a total 
of 14 periodic orbits up to period 3 along the same proce- 
dure adopted in [^7j. Resonances belonging to different 
symmetry characters can be identified in the full-space 
configuration and compared with the semiclassical pre- 
dictions J27]]. A similar comparison can be performed 
in the fundamental domain, only based on resonances in 
the B2 subspace. As for the 3-disk case, the semiclassi- 
cal theory provides a qualitatively fair prediction of the 
resonance frequencies. We observe that, in general, the 
agreement for the lowest-energy widths is not as satisfac- 
tory, with discrepancies increasing with decreasing fre- 
quencies. This kind of discrepancy, which is also found 
in the 3-disk geometry discussed above, is intrinsic to the 
semiclassical calculation because of the large correction of 
the stationary phase approximation [ffijf . For the 3-disk 
|p|, p^0|j3l| and 4-disk J27j systems, where exact quantum- 
mechanical calculations are available, the very low-lying 
resonance widths of semiclassical resonances appear to 
be systematically smaller (up to a factor 3) compared to 
the corresponding quantum ones. 

A few general remarks are in order. Although the 
agreement between the experimental scattering reso- 
nances and the corresponding semiclassical predictions 
is generally within a few percents (sa 5%), some dis- 
crepancies are also evidenced from the data we analyzed. 
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Such discrepancies may manifest in the form of both fre- 
quency shifts or width modifications of the predicted res- 
onances as well as in the presence of additional peaks in 
the experimental trace. Various mechanisms and experi- 
mental limitations are expected to contribute as possible 
sources of errors, including: symmetry-breaking pertur- 
bations introduced by non-perfect geometries, effects as- 
sociated with spurious reflections, non-idealities in the 
operations of the microwave absorbers (e.g., frequency- 
dependent response), or slight height variations over the 
cavity area. The combined action of such mechanisms 
makes open-geometry microwave experiments compara- 
tively more demanding with respect to their closed-cavity 
counterpart, where some of the above error sources are 
practically irrelevant. While a deeper understanding of 
the unavoidable non-idealities faced by the experiments, 
along with the necessary technical improvements, are 
likely to be necessary for establishing a fully quantita- 
tive detailed comparison, the level of agreement reached 
in our present investigation can be considered a very sat- 
isfactory match with the opportunity of retaining a rela- 
tively simple experimental methodology. 



center point ko has been checked in the calculations. We 
consider the average of autocorrelations with a different 
fco to compensate for slight dependences on the center 
point. Plots of typical experimental autocorrelations for 
2-, 3-, and 4-disk systems are shown in Fig. B, Fig. ^ 
and Fig. ^ respectively. 

By inserting the explicit representation of |5*2i(fc)| 2 as 
a sum of Lorentzians, Eq. (|30|), the autocorrelation is 
found as 



CiCjiji+Jj) 



(k - (s l 



Sj)) 2 + hi + 1 3 ) 2 



(33) 



In a regime where there are no overlapping resonances, 
\ s i ~ s j\ ili + 7j); an d t ne small- k behavior of the 
autocorrelation can be simplified as [pi] 



C(k) 



c 2 + 4 7 2 



(34) 



By exploiting a result from semiclassical Random Ma- 
trix Theory, the above sum can be replaced by a single 
Lorentzian |37|J38il: 



VI. SPECTRAL AUTOCORRELATION 

We now turn to analyze the data in terms of the so- 
called spectral autocorrelation function, which was cal- 
culated as 

C( K ) = {\S 21 (k-( K /2))f\S 21 (k+(n/2))f) k . (31) 

Here, k is the wave- vector difference and ( ) k denotes 
an average over a band of wave- vectors centered at some 
value k = ko and of width Ak, the latter being large 
enough to include an appreciable number of resonances. 
The average (||l]) also includes a suitable window function 
which is chosen as |35| 



C(k) = C(0) 



1-N/V6 







Id > V6 



k — kg 
Ak 



(32) 



Besides its intrinsic interest, an additional motivation 
for investigating the properties of C(n) comes from the 
correspondence, pointed out in Sec. IIIc, with experi- 
ments performed on mesoscopic transport. For ballis- 
tic conductors, a formally similar magnetic-field correla- 
tion function received extensive theoretical and experi- 
mental attention as a potential probe for quantum chaos 
jl0|,[35]]. A similar autocorrelation measure was also con- 
sidered recently in the context of molecular photodisso- 
ciation spectra p6| . In our microwave experiments, the 
wave- vector plays the role of the magnetic field and since 
|5 , 2i(^)| 2 oc T{k), the function C(k) can be regarded as 
a measure of the wave-vector correlations of the two- 
probe conductance. The dependence of the autocorre- 
lation function on both the finite window Afc and the 



1 + [Khf 



(35) 



where the parameter 7 = 7^ is identified with the clas- 
sical escape rate from the chaotic scattering region, with 
the velocity scaled to 1. Accordingly, one can interpret 
the width of the autocorrelation function as an average 
width (and thereby lifetime) of the resonances |59],flO| . 



A. Universal features: classical escape rate 

According to the above predictions, a universal behav- 
ior of the the autocorrelation function is expected for 
sufficiently small correlation scales, regardless of the de- 
tails of the geometry and the way the system is excited. 
Such universal behavior is captured by the single classical 
parameter 7. We recall its definition. Classically, if we 
shoot particles towards the scatterer, the number N(t) 
of particles remaining in the scattering region after time 
t decays exponentially as 



N(t) = N(0) exp(-7cj t) . 



(36) 



where %i = A(l — d) is the classical escape rate, A is the 
Lyapunov exponent of the manifold of infinitely trapped 
orbits (strange repeller), and d is the information dimen- 
sion of the unstable manifolds. The scaled escape rate, 
corresponding to unit velocity, is defined as j c i = jd/v, v 
being the speed of the particles. The classical escape rate 
can be calculated through the classical Ruelle ^-function 

11, 
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C/3(s)=n [l-(-l) L "exp(* 



(37) 



which is analytical in the half-plane Res < —P(j3), and 
has poles in the other half-plane. In particular, £p(s) has 
a simple pole at s = P{f3). Here, is the so-called 

Ruelle topological pressure, from which all the charac- 
teristic quantities of classical dynamics can be derived in 
principle. The classical escape rate is 7c ; — —P(l). 

For the various scattering geometries investigated ex- 
perimentally, we calculated the appropriate autocorre- 
lation function from the observed trace and fitted the 
small- « portion of the resulting curve with the Lorentzian 
behavior ( |35|) , thereby extracting an experimental escape 
rate 7 ?m . In general, good agreement is observed with the 
classical prediction 7c /, implying that in the regime of 
universality the characteristic scale of wave-vector corre- 
lations in the measured two-point quantum conductance 
is well reproduced by knowledge of the chaotic classical 
scattering dynamics, through the classical escape rate *y c i . 

2- disk: For the integrable two-disk system, the infor- 
mation dimension d = 0, thus j c i = A. For unit velocity 
and R > 2a, we get pi 



7c/ = A = 



1 



R-2a 



In A 



(38) 



A = A(er) being the eigenvalue of the monodromy ma- 
trix introduced in Eq. J24|). The autocorrelation for the 
experimental set up with a — 5 cm, R = 40 cm (a = 8) 
is shown in Fig. 5. A value j qm = 0.083 cm -1 is found, 
which is in excellent agreement with the classical result 
7c/ = 0.088 cm" 1 . 

3-disk: For the three-disk scatterer, we have flq] 



7c/ ~ 4 In (1.072 a; 



(39) 



A representative wave- vector autocorrelation function for 
this system is displayed in Fig. 3, where the funda- 
mental domain configuration has been investigated with 
a = 5 cm and R = 20 V3 cm. From equation (39), 
7c/ ~ 0.058 cm -1 , to be compared with the experimental 
value 7om = 0.064 cm -1 . The latter, leading to a scaled 
value 7 gm a = 0.32, is in very good agreement with both 
semiclassical and Monte Carlo estimates as given by . 
The quantitative agreement between the small- k decay of 
correlations and the Lorentzian curve demonstrates ex- 
plicitly that behavior in such a region is universal, with 
no dependence on the actual details of the geometry. 

4-disk: For the four-disk scatterer, we use the results 
of Refs. for comparison with the experimental data. 
Asymptotically, for large R, 



1 



7c/ 



V2R 



ln(2\/2cr) 



(40) 



The autocorrelation function for 4-disk data in the fun- 
damental domain is displayed in Fig. 7, corresponding 



to a = 5 cm and R — 20\/2 cm. The value 7gm = 0.070 
cm -1 compares pretty well with the estimate from Eq. 
©, 7c/ = 0.069 cm" 1 . 

In Fig. 8, the experimental escape rates 7 qm of the 
4-disk system are compared with the classical prediction 
7c/ for several values of the ratio a. Note that data are 
included for 17 configurations of the different reduced 
(1/8, 1/4, 1/2 and full space) representations of the 4- 
disk geometry shown in Fig. 4. The radius of the disks 
used was a — 5 cm for the 1 / 8th space, and a = 2 cm for 
the others. A relevant quantity is the abscissa of absolute 
convergence s c for Eq. (|l7]) , which can also be estimated 
from the Ruelle (^-function with the classical cycle weights 
t p replaced by the corresponding semiclassical ones. s c 
serves as a crude lower bound of the escape rate ^] . The 
latter is also shown in Fig. 8. s c becomes negative for 
R/a < 4.5. 

By comparing the results found for the escape rate 
while passing from 2- up to 4-disk scattering, progres- 
sively smaller values are obtained. In general, it is inter- 
esting to examine the variation of the escape parameter 
with increasing number n of scatterers. For n — > oo, one 
obtains a so-called Lorentz scatterer ||. We carried out 
experiments with n = 20, with corresponding (scaled) 
escape rate 7gm ~ 0.05, which is roughly an order of 
magnitude smaller than the three- or four-disk values. 
In agreement with physical intuition, this indicates that 
the system approaches a closed system when the number 
of disks becomes very large. Accordingly, the escape rate 
from the chaotic region is found to be quite small. 



B. Non-universal features 

For intermediate k, the semiclassical prediction of Eq. 
(|35| ) fails because of the presence of the periodic orbits, 
which leads to non-universal behavior. In the case of just 
one periodic orbit, one may express the full two-point 
correlation function as 



C(k) cx 



27 



^ (k- 

n=0 v 



nAs) 2 + 47 2 



(41) 



where As denotes the spacing between resonances in the 
wave- vector domain. For example, for the two-disk prob- 
lem discussed above (Sec. IVa), As = 2ir/(R— 2a), 
provided the space is probed through the antisymmetric 
^-representation only (32|. Thus, the autocorrelation 
oscillates with period As. Very good agreement is found 
between experiment and theory for this integrable two- 
disk system (Fig. 5), where the expected value As ~ 0.21 
cm -1 is identical with the observed oscillation period. 

For the three-disk scatterer in the fundamental do- 
main, the average length of the periodic orbit per period 
is roughly given by (Iq + h)/2, the mean separation of 
the resonances being therefore 
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Asfd = 



-Itt 



2R-(2 + V3)a 



(42) 



The autocorrelation is expected to oscillate with a period 
roughly equal to Asfd- The value Asfd — 0.25 cm -1 
predicted from Eq. (42) is in very good agreement with 
the scale of the oscillations in Fig. 6. 

Let us finally discuss the four-disk scatterer. In the full 
space, the average length of the periodic orbits per pe- 
riod can be estimated as the average length of the eight 
periodic orbits, 12, 23, 34, 41, 1234, 1432, 13, 24, where 
1,2,3,4 are the labels of the four disks |2(J. The mean 
separation between the resonances is then given approx- 
imately by 



As 



2tt 



full 



(2 + l/V2)R - (3 + V2)c 



(43) 



For the 4-disk system in the fundamental domain (one- 
eighth of the phase space), the average length of the pe- 
riodic orbit per period is (Iq + ^)/3, thus the mean 
separation is 



As 



Gtt 



FD 



(3 + V2)R - 2(2 + V2)c 



(44) 



The autocorrelation will oscillate with an approximate 
period equal to Asfd, which indicates the deviation 
from the semiclassical theory due to the presence of 
the periodic orbits. Thus, the large wave-vector (or 
short time) behavior is system specific. The value of 
Asfd = 0.21 cm -1 found from Eq. ( |4^ ) is in good agree- 
ment with the scale of oscillations in Fig. 7. 

Non-universal contributions can play in general a cru- 
cial role in determining the overall structure of the spec- 
tral autocorrelation, since they can be of the same or- 
der of the universal result of Random Matrix Theory. 
Semiclassical methods have recently provided an insight- 
ful tool in modeling non-universal properties in addi- 
tion to universal ones Beside the general remarks 
mentioned above, the systematic identification of non- 
universal features, along with their interplay with the 
universal scattering properties, deserves a separate ex- 
perimental investigation. 



VII. DISCUSSION AND CONCLUSIONS 



Two main conclusions can be drawn from the exper- 
iments. First, the general validity and the predictive 
power of semiclassical methods have been tested directly 
by comparing the observed spectra with the correspond- 
ing semiclassical predictions. Qualitative agreement has 
been verified in all the situations investigated, and quan- 
titative comparison found for a wide class of relatively 
sharp resonances. In addition, the experiments point 
out the validity of the wave-vector autocorrelation func- 
tion as a probe for phase-space structure and quantum 
chaos. Values of wave-mechanical escape rates have been 
extracted from the observed autocorrelations, and com- 
pared to their classical counterparts. In general, the 
agreement observed between the measured escape pa- 
rameters and the corresponding semiclassical predictions 
tends to be quantitatively superior than the one reach- 
able in the detailed comparison of single resonance peaks. 
A similar conclusion has been reached for closed geome- 
tries where a much better agreement between experimen- 
tal and numerical results has been obtained for the sta- 
tistical properties rather than for the comparison of the 
individual resonances |^]. This feature reflects the na- 
ture of the escape rate as an average spectral property. 

We remark that the present experiments, which 
probe wave-vector dependence, nicely complement mea- 
surements performed on semiconductor microstructures, 
where a similar role has been stressed for the magnetic 
field correlation of the two-point conductance. From the 
broader perspective of quantum-classical correspondence, 
the experiments shed light on the interplay between clas- 
sical and quantum features of the scattering dynamics, by 
showing that measurable properties like quantum corre- 
lation lengths can be predicted from a knowledge of the 
classical chaotic scattering behavior. 

Our investigation clearly points out, among other is- 
sues, the need for a deeper understanding of the non- 
universal properties of the spectral statistics and their 
interplay with universal ones. Such an investigation is 
likely to involve wave- vector autocorrelation functions as 
considered in the present analysis, as well as different 
quantities useful to characterize and probe the chaotic 
dynamics. In general, novel tools may be required in or- 
der to pull out the whole amount of information encap- 
sulated in the observed spectra. Work is ongoing along 
these directions. 



We presented an extensive experimental investigation 
of hard-disk chaotic scattering in microwave open cavi- 
ties. The experiments provide a conceptually clean and 
direct realization of the n-disk open billiard problem. By 
exploiting our (in principle unlimited) capability to vary 
the geometry, chaotic scattering was studied in various 
configurations by both changing the number of disks and 
the symmetry properties of the underlying phase space. 
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FIG. 1. Photo of the experimental apparatus. (Top) 
Closed cavity. (Bottom) Open cavity. 
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FIG. 2. Experimental (solid line) vs. semiclassical 
(dashed line) transmission function 15*21 1 2 for a 2-disk system 
with R = 40 cm and a — 5 cm probed in the fundamen- 
tal domain. (Inset) Sketch of the corresponding experimental 
configuration. The separation distance R is indicated. 
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Experimental (solid line) vs. 



20 



semiclassical 



(dashed line) transmission function IS21I for a 3-disk system 
in the fundamental domain with R = 20\/3 cm and a — 5 cm. 
The corresponding experimental configuration is sketched in 
the inset. The separation distance 7? is indicated. 
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FIG. 4. Experimental (solid line) vs. semiclassical (dashed 
line) transmission function |S2i| 2 for a 4-disk system. (Top) 
Full space geometry with R = 8 cm and a = 2 cm. (Bot- 
tom) l/8th space (fundamental domain) with R = 20 cm and 
a — 5 cm. The corresponding experimental configurations are 
sketched in the insets. The separation distance R is indicated. 
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FIG. 5. Wave- vector autocorrelation C{k) of the 2-disk sys- 
tem with R=40 cm and a=5 cm. Data are shown for the 
half-space configuration of the 2-disk geometry, correspond- 
ing to the A2 representation. The correlation is calculated 
with interval Afe = 3 cm -1 . The different sets represent dif- 
ferent values of the central wave- vector ko- The bold line is a 
Lorentzian with 7 9m =0.083 cm -1 . 
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FIG. 6. Wave- vector autocorrelation C(k) of the 3-disk sys- 
tem with R — 20\/3 cm and a = 5 cm. Data are taken in 
the fundamental domain, corresponding to the A2 represen- 
tation. The correlation is calculated with interval Afc = 2 
cm -1 . The different sets represent different values of the 
central wave- vector fco- The bold line is a Lorentzian with 
7 gm = 0.064 cm" 1 . 
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FIG. 7. Wave- vector autocorrelation C(k) of the 4-disk sys- 
tem with R — 20\/2 cm and a = 5 cm. Data are taken in 
the fundamental domain, corresponding to the B2 represen- 
tation. The correlation is calculated with interval Afc = 3 
cm -1 . The different sets represent different values of the 
central wave- vector fco. The bold line is a Lorentzian with 
7 qm = 0.070 cm" 1 . 
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FIG. 8. Experimental escape rate 7 9m scaled to radius 
a = 1 vs. ratio a = R/a. Data are shown for different 
reduced configurations of the 4-disk geometry: l/8th space 
(open squares), 1/2 space (open circles), l/4th space (filled 
squares), full space (triangles). The classical escape rate (solid 
line) is calculated from the first 3 periodic orbits in the funda- 
mental domain. The abscissa of convergence s c of Eq. (17), 
which is shown as a dot-dashed line, represents a lower bound 
on the quantum escape rate. 
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